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Abstract The relation between the mass of supermassive black holes located in the 
center of the host galaxies and the kinetic energy of random motions of the correspond- 
ing bulges can be reinterpreted as an age-temperature diagram for galaxies. This rela- 
tion fits the experimental data better than the Mm — Mq, Mm — Lq, and Mm — a laws. 
The validity of this statement has been confirmed by using three samples extracted 
from different catalogues of galaxies. In the framework of the /1CDM cosmology our 
relation has been compared with the predictions of two galaxy formation models based 
on the Millennium Simulation. 
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1 Introduction 

Nowadays an increasingly number of multi-band observational evidences reinforces the 
picture that all massive galaxies, with a spheroidal component, contain a super massive 
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black hole (SMBH, Mm > 10 6 Mq) at their center. Some kind of connection between 
the traits of central SMBHs and the evolutionary state of their host galaxies should 
be expected and actually it is currently under study and debate. In particular, the 
mass of the SMBHs seems to be closely related with the properties of the spheroids 
in which they reside. Many efforts have been made in this sense in order to achieve a 
pragmatic correlation, such as for example the ones between the mass of the SMBH and 
the bulge 1 stellar mass or luminosity (Mm — Mq; Mm — Lq) [1-7], velocity dispersion 
(Mm - a) [8-10], effective radius (Mm - Re) [6], kinetic energy (Mm - M G a 2 ) [11,12], 
Sersic index (M, — n) [13,14]. However, as already remarked by Novak et al. [15], it 
is difficult to understand the fundamental nature of the correlations between SMBH 
and host properties, because all such relations depend critically on the accuracy of the 
published error estimates in all quantities under consideration. Moreover, the conser- 
vative practice of adopting large error bars turned out not to be such a good strategy 
if only few data are considered for the fit. 

A more complex approach is to read the above-mentioned relations in terms of a 
fundamental plane, which relates the SMBH mass to two or more spheroid properties 
such as galaxy effective radius, stellar velocity dispersion, luminosity, and mass [9,16]. 
A hydrodynamical simulation of major galaxy mergers, including the effects of black 
hole (BH) accretion and feedback, supports this picture [17]. Other authors reached 
the same conclusion by modeling the cosmological co-evolution of galaxies and their 
central SMBHs within a semianalytical framework [18]. 

The possibility that the mass of a SMBH correlates with the total gravitational 
mass of its host galaxy, or with the mass of the dark matter halo in which it pre- 
sumably formed, has been investigated [19,20], and supported by recent self-consistent 
simulations of the co-evolution of the SMBH and galaxy populations [21], which asserts 
that the mass of a SMBH is determined primarily by the mass of the dark matter halo. 

Also important is to understand the mechanisms that regulate the growth rate of 
the SMBH. Studies of the distribution of Eddington luminosity ratios, ibol/^Edd) of 
active galactic nuclei (AGNs) show that the energy-storing rate in luminous AGNs are 
ultimately determined by SMBH self-regulation of the accretion flow and this process 
acts on large scales (> 1 kpc) [22]. 

An alternative fundamental plane has been proposed in order to correlate the radio 
and the X-ray luminosity, and the BH mass [23,24]. However, this plane turns into an 
effective BH mass predictor, like the very popular M, — a relation, only if the obscured 
AGNs are excluded [25]. Actually, it is not clear at all if the intrinsic activity (Seyfert, 
Liner, etc.) of galaxies gives them a special place in the fundamental planes as well as in 
any diagram originated by the other proposed relations. This topic definitely deserves 
a deep investigation. 

Very recently, Feoli & Mancini [26] (thereafter FM) suggested a fascinating con- 
nection between the stored energy of central SMBHs and the evolutionary process of 
galaxies. Using a sample of 64 galaxies, they investigated the relation between the mass 
of the SMBHs and the kinetic energy of the random motion of the corresponding galaxy 
bulges. Besides the fact that, as already noted in previous papers [11,12], this relation 
works better than the most common Mm — a and M. — Mq laws, they found some 
analogies between the Mm — Mqo 2 /c 2 plane for galaxies and the Hertzsprung-Russell 
(HR) diagram for stars. The HR diagram connects the energy radiated by the nucleus 



1 Here bulge refers to either the spheroidal component of a spiral/lenticular galaxy or to a 
full elliptical galaxy. 
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of a star with its surface temperature. Similarly, the FM diagram essentially connects 
a property of the inner nucleus of a galaxy, the energy stored by the SMBH, M»c 2 , 
with a property of the external surface of its bulge, i.e., the kinetic energy of random 
motions. Moreover, each morphological type of galaxy generally occupies a different 
area in the FM plane, reminiscent of the different positions occupied by stars of the 
various spectral classes in the HR diagram. 

In this paper, we want to verify the goodness of the M. — Mqv 2 relation, first 
of all comparing it with the other popular relations (essentially M. — Mq, M. — Lq, 
and M, — a). In order to do this, we consider three samples of galaxies extracted from 
three different catalogues. Thanks to this procedure, we avoid the possibility that the 
tightness of the FM relation might depend on a suitable choice of the data. Then, we 
complete our analysis comparing the results, derived by real data, with the predictions 
of two semi-analytic hierarchical models that follow the cosmological co-evolution of 
dark matter, galaxies and black holes. 

The paper is structured as follows: in § 2 we provide our reference catalogues and 
report the data of the correspondent samples. In § 3 we discuss the results which 
emerge from our analysis. These results are then compared with the predictions of 
two hierarchical models of galaxy structure formation in § 4. Finally, § 5 contains the 
summary of the main conclusions of this work. 

2 The samples 

In order to check which relation is the most effective to predict the black hole mass, 
we have tested the M, — a, the M, — Mq, the M» — Lq and the M, — Mqo 1 relations 
on three different samples of galaxies extracted from the catalogues of Graham [27], 
Giiltekin et al. [28], and Hu [29] respectively. Since these three catalogues have been 
compiled independently, we are quite safe from any tampering of the data due to a 
personal intervention, that amounted to selection from the literature of the best values 
for the mass of the galaxies, which are missing in the first two catalogues. 
The relations that we want to study can be written in the following form 

lo Slo ( M ») = b + mlog w x, (1) 

where m is the slope, b is the normalization, and a; is a parameter of the host bulge 
such as the mass (Mq), the luminosity (Lq), the central velocity dispersion (a c ), 
the effective velocity dispersion (<r e ), or the kinetic energy of the random motions 
(Mqcj 2 / (?). Eq. 1 can be used to predict the values of M. in other galaxies once we 
know the value of x. In order to minimize the scatter in the quantity to be predicted, 
we have to perform an ordinary least-squares regression of M. on x for the considered 
galaxies of which we already know both the quantities. In Table 1 we collect the fits 
obtained for all the relations for each sample. As in [12,14,26], these fits were obtained 
taking into account the error bars in both variables and using the routine FITEXY [30] 
for the relation y = b + mx, by minimizing the \ 2 ■ 

2.1 Graham sample 

Graham (2008) compiled a catalogue of 76 galaxies with direct SMBH mass measure- 
ments [27]. We select only 61 galaxies (27 ellipticals, 19 lenticulars, 15 spirals), since 
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a reliable value for the central velocity dispersion and the mass was not available for 
all the host bulges of the catalogue. These galaxies are listed in Table 2 together with 
the values of their respective parameters. Concerning the errors in the measures, we 
adopt the same strategy as in [11,12]: We consider that the error for the bulge mass 
is 0.18 dex in log 10 MG for all the galaxies [7], while the relative error on the velocity 
dispersions is 5%. 

In Fig. l(a)-(c), we reported the M, — a, M, — Mq, and M, — Mqo 2 relations in 
log-log plots (we associated a particular marker to each galaxy according to its morpho- 
logical type). The best-fitting lines are also shown for each diagram. The comparison 
of the fits of the three relations reveals that the x 2 > the intrinsic scatter so, and the 
Pearson linear correlation coefficient r of the M» — Mqo 2 relationship are better than 
the other ones (Table 1). 

The lenticular galaxies NGC4342 and NGC4350 have been excluded from the fits 
because the SMBH mass to bulge mass ratio for these galaxies is very high [37,51]; 
their positions in Fig. 1(c) are enclosed in a dashed box. In Fig. 1(c) we also marked the 
position of the Milky Way and of the spiral NGC4594 (the Sombrero galaxy), which is 
surrounded by a halo of stars, dust, and gas that indicate it may actually be described 
as an elliptical galaxy that contains a more robust interior configuration [50]. We refer 
the reader to [26] for a comprehensive description of the diagrams reported in Fig. 1. 

2.2 Giiltekin et al. sample 

In their catalog, Giiltekin et al. (2009) listed 55 galaxies, that had a direct measure- 
ment of the mass of their SMBHs, and gave the upper limits to SMBH masses of 18 
additional galaxies; the latter are not considered in our analysis [28]. The values of the 
parameters of the 55 galaxies (27 ellipticals, 14 lenticulars, and 14 spirals) are listed in 
Table 3. Again we consider that the error for the bulge mass is 0.18 dex in log 10 MG, 
whereas we use an absolute error of ±0.10 mag for the visual magnitudes there where 
it is absent for some galaxies of the catalogue. We consider only 52 galaxies in the 
study of the M. — Ly relation, since the values of the magnitude of three of them 
have not been reported by Giiltekin et al. in their catalogue [28]. Moreover, these au- 
thors preferred to report the values of the effective dispersion velocity instead of the 
central one. However, as already noted by different authors [15,56], the two ways of 
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Table 2 The data are taken from Graham (2008) [27], except for the values of Mq whose 
references are given in the last column. The value of cr c for Abcll 1836 is taken from [28] . The 
error for the bulge mass is 0.18 dex in the log 10 Mq for all galaxies, while the relative error 
on the velocity dispersions is 5%. Brightest cluster galaxy (BCG) is defined as the brightest 
galaxy in a cluster of galaxies 
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0.7 


x 


10 7 


1.1 x 10 11 


[47] 


NGC4258 


SB 


134 


3.9 


X 


10 7 
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1 2 3 4 5 6 7 



Log |0 [M G cr 2 /M G c 2 ] 

Elliptical galaxies Dwarf elliptical galaxies 

^ Lenticular galaxies £ Barred lenticular galaxies 

Spiral galaxy Barred spiral galaxies 

Fig. 1 Best-fitting (a) M m — <x, (b) M. — Mq, and (c) M. — Mqcj 1 relations for the sample 
extracted from Graham (2008) [27]. The symbols represent elliptical galaxies (red ellipses), 
lenticular galaxies (green circles), barred lenticular galaxies (dark green circles), spiral galaxy 
(blue spirals), barred spiral galaxies (dark blue barred spirals), and dwarf elliptical galaxies 
(orange ellipses). The dashed box encloses the lenticular galaxies NGC4342 and NGC4350, 
which have been excluded from the fits due to their high value of the M, /Mq ratio. 
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Table 3 The data are taken from Gultekin et al. (2009) [28], except for the values of Mq 
whose references are given in the last column. The error for the bulge mass is 0.18 dex in the 
log 10 Mq for all galaxies. 



Galaxy 


Type 


a c (km/s) 




M v 






M, 


(Mq) 


SM, 


• (Mq) 


Mq(Mq) 


References 


A1836-BCG 


E 


288 ± 14 


—23. 


31 


± 


0. 
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[29] 


A3565-BCG 


E 
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.8 
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E 
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0. 


10 


2 


,7 


X 
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,7 


X 
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[6] 
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E4 
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± 


0. 


15 


2 


.8 
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.2 
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[6] 
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75 ± 3 
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83 
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0. 


05 


3 


1 
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6 


X 
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[7] 


NGC821 


E4 
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—21. 


.24 


± 


0. 


.13 


4 


.2 


X 
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3 


.5 


X 
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1.3 x 10 11 


[7] 


NGC1399 


El 
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—22 


13 


± 





10 


1 


3 


X 


10 9 


6 


6 
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[32] 
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[33] 
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[7] 
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.1 


X 


10 8 


9.7 x 10 10 


[7] 


NGC4261 


E2 


315 ± 15 


—22 


72 


± 


0. 


06 


5 


.5 


X 


10 s 


1 


2 


X 


10 s 


3.6 x 10 11 


[7] 


NGC4291 


E2 


242 ± 12 


-20. 


67 







13 


3 


2 


y 


10 8 


3 


x 


x 


10 8 


1.3 x 10 11 


[7] 


NGC4374 


El 


296 ± 14 


-22 


45 







.05 


x 


.5 


V 


10 9 


x 


x 


x 


10 9 


3.6 x 10 11 


[7] 


NGC4459 


E2 


167 ±8 


-21. 


06 




0. 


.04 


7 


4 


V 


10 7 


x 


4 


x 


10 7 


7.9 x 10 10 


[33] 


NGC4473 


E4 


190 ±9 


-21. 


14 




0. 


.04 




3 


x 


10 8 


o 


.9 


x 


10 8 


9.2 x 10 10 


[7] 


NGC4486 


El 


375 ± 18 


-22 


92 







.04 


3 


.6 


x 


10 9 


x 





x 


10 9 


6.0 x 10 11 


[7] 


NGC4486A 


E2 


111 ± 5 


-18. 


70 







05 


x 


3 


X 


10 7 





.5 


x 


10 7 


4.1 x 10 9 


[34] 


NGC4649 


E2 


385 ± 19 


-22. 


.65 







.05 


2 


x 


X 


10 9 


o 


.6 


x 


10 9 


4.9 x 10 11 


[7] 


NGC4697 


E6 


177 ±8 


-21 


29 


± 


0. 


.11 


2 


.0 


x 


10 8 





2 


x 


10 8 


1.1 x 10 11 


[7] 


NGC4742 


E4 


90 ± 5 


-19. 


.91 


± 


0. 


10 




.5 


X 


10 7 





6 


x 


10 7 


6.2 x 10 9 


[7] 


NGC5077 


E3 


222 ± 11 


-22. 


04 




0. 


13 


g 


o 




10 s 


5 


o 




10 8 


2.1 x 10 11 


[35] 


NGC5576 


E3 


183 ±9 


-21. 


.26 


± 





13 


X 


.8 


V 


10 8 





x 


x 


10 8 


1.5 x 10 11 


[54] 


NGC5845 


E3 


234 ±11 


-19. 


77 


± 





13 


2 


.9 


y 


10 s 


x 


7 


x 


10 s 


3.7 x 10 10 


[7] 


NGC6251 


El 


290 ± 14 


-22 


90 


± 


0. 


10 


6 


.0 


V 


10 8 


2 


.0 


x 


10 8 


5.6 x 10 11 


[7] 


NGC7052 


E3 


266 ± 13 


-22. 


.35 


± 





10 


4. 


.0 


X 


10 8 


2. 


.8 


X 


10 8 


2.9 x 10 11 


[7] 


NGC3115 


SO 


230 ± 11 


-21. 


18 




0. 


.05 


9 


.6 




10 s 


5 


4 




10 8 


1.2 x 10 11 


[7] 


NGC3245 


SO 


205 ± 10 


-20. 


96 


± 


0. 


10 


2 


2 


x 


10 s 





.5 


x 


10 8 


6.8 x 10 10 


[7] 


NGC3585 


SO 


213 ± 10 


-21. 


.80 


± 


0. 


20 


3 


4 


x 


10 s 


x 


.5 


x 


10 8 


1.8 x 10 11 


[54] 


NGC3998 


SO 


305 ± 15 


-20. 


32 


± 


0. 


10 


2 


x 


x 


10 8 


2 


x 


x 


10 8 


5.5 x 10 10 


[36] 


NGC4026 


SO 


180 ±9 


-19. 


83 




0. 


20 


2 


x 




10 8 


o 


7 




10 8 


5.2 x 10 10 


[7] 


NGC4342 


SO 


225 ± 11 


-18. 


84 




0. 


10 


3 


Q 




10 s 


2 


Q 




10 8 


1.2 x 10 10 


[7] 


NGC4564 


SO 


162 ±8 


-19. 


60 




0. 


32 


(5 


.9 


x 


10 7 


x 


.0 


x 


10 7 


4.4 x 10 10 


[7] 


NGC5128 


SO/E 


150 ±7 


-21 


82 




0. 


08 


3 







10 8 


o 


4 




10 8 


3.6 x 10 10 


[38] 


NGC5252 


SO 


190 ± 10 












X 


o 




10 9 


x 


.6 




10 9 


2.4 x 10 11 


[6] 


NGC7457 


SO 


67 ±3 


-18 


72 


± 


0. 


11 


x 


x 


x 


10 6 


x 


7 


x 


10 6 


7.0 x 10 9 


[7] 


NGC1023 


SBO 


205 ± 10 


-20. 


.61 


± 


0. 


.28 


x 


.6 


x 


10 7 





.5 


x 


10 7 


6.9 x 10 10 


[7] 


NGC2787 


SBO 


189 ±9 


-18 


90 







10 


4 


.3 




10 7 


Q 


5 




10 7 


2.9 x 10 10 


[40] 


NGC3384 


SBO 


143 ± 7 


-19. 


93 







22 


x 


g 




10 7 


Q 


.3 




10 7 


2.0 x 10 10 


[7] 


NGC4596 


SBO 


136 ±6 


-20 


70 


± 





.10 


8. 


A 


X 


10 7 


4. 


.4 


X 


10 7 


2.6 x 10 10 


[6] 


Oircinus 


c 
o 


1 f^fi _|_ 1 Q 

loo it lo 


"1 1 


36 


± 


u. 


i n 


1 


.7 


X 


10 6 





.4 


X 


10 6 


o.U X 1U 


\A 1 1 

[41j 


Milky Way 


s 


105 ± 20 












4 


1 


X 


10 6 





.6 


X 


10 6 


1.1 x 10 10 


[7] 


NGC224 


s 


160 ±8 


-21 


84 


± 





30 


1 


.5 


X 


10 8 





.9 


X 


10 8 


4.4 x 10 10 


[42] 


NGC1068 


s 


151 ± 7 


-22 


17 


± 


0. 


10 


8. 


6 


X 


10 6 





.3 


X 


10 6 


1.5 x 10 10 


[43] 


NGC2748 


s 


115 ± 5 


-20. 


97 


± 


0. 


10 


4. 


,7 


X 


10 7 


3 


.8 


X 


10 7 


1.7 x 10 10 


[44] 


NGC3031 


s 


143 ±7 


-21. 


.51 


± 


0. 


10 


8. 


.0 


X 


10 7 


6 


.9 


X 


10 7 


1.0 x 10 10 


[45] 


NGC4594 


s 


240 ± 12 


-22. 


44 


± 


0. 


15 


5 


.7 


X 


10 s 


5 


.3 


X 


10 8 


2.7 x 10 11 


[7] 


NGC4945 


s 


134 ±7 












1 


.4 


X 


10 6 





,7 


X 


10 6 


3.0 x 10 9 


[41] 


NGC1300 


SB 


218 ± 10 


-21. 


34 


± 


0. 


10 


7. 


.1 


X 


10 7 


3 


.6 


X 


10 7 


2.1 x 10 10 


[44] 


NGC3227 


SB 


133 ± 12 


-20. 


73 


± 


0. 


10 


1 


.5 


X 


10 7 





.8 


X 


10 7 


3.0 x 10 9 


[5] 


NGC4151 


SAB 


93 ±5 


-20. 


68 


± 


0. 


10 


4. 


.5 


X 


10 7 





.5 


X 


10 7 


1.1 x 10 11 


[47] 


NGC4258 


SAB 


115 ± 10 


-21. 


.31 


± 





10 


3 


.8 


X 


10 7 





.1 


X 


10 6 


1.1 x 10 10 


[6] 


NGC4303 


SAB 


84 ±4 


-21. 


65 


± 


0. 


10 


4. 


.5 


X 


10 6 


9 


.5 


X 


10 6 


1.6 x 10 9 


[55] 


NGC7582 


SB 


156 ± 19 


-21. 


.51 


± 


0. 


10 


5 


.5 


X 


10 7 


1 


.6 


X 


10 7 


1.3 x 10 11 


[49] 
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(a) (b) 




Log lo [Z.y/L ] Log I0 [M G cr 2 /M G c 2 ] 

Fig. 2 Best-fitting (a) M. - a, (b) M. - M G , (c) A/. - L v , and (d) Af. - M G a 2 relations 
for the sample extracted from Gultekin et al. (2009) [28]. The symbols are the same of Fig. 1. 

measuring the velocity dispersion do not generate profound differences. The V-band 
luminosities have been calculated from the extinction-corrected magnitudes, -My ^ u i ge , 
using log 10 (l/v/i0,v) — 0.4 (4.83 — M Vi bulgc)i m accordance with [28]. 

The diagrams of the M» — a, M» — Mq, Aim — Ly, and M» — Mqo 2 relations are 
shown in log-log plots together with the best-fitting lines, see Fig. 2(a)-(d) (the symbols 
are the same of Fig. 1). From inspection of Table 1, it is clear that the M. — Mqo 1 
relation works better than the other relations when applied on the same sample of 
galaxies. 



2.3 Hu sample 

Hu (2009) compiled a catalogue of 58 galaxies (28 ellipticals, 19 lenticulars, and 11 spi- 
rals) [29]. The main parameters are summarized in Table 4. He estimated the dynamical 
mass of bulges by Mq ^ yn = kR e a c /G, where Re is the bulge effective radius, G is 
the gravitational constant, and k is a function of the Sersic index which is determined 
numerically. 



9 




Of course, there are also other methods to estimate the mass of the galaxies, and in 
particular we prefer those which use the Jeans equation or three-dimensional models, 
see appendix B of [26]. 

The la errors of Lk and °"c are adopted as 10% and 5% respectively, whereas that 
of logM G is adopted as 0.15 dex (40%) [29]. 

In Fig. 3(a)-(d), we show the Mm — a, Mm — M G , M, - L K , and M, - M G a 2 
relations in log-log plots (the symbols are the same as those used in Fig. 1), together 
with the best-fitting lines. We recovered the same result of [29], that is the intrinsic 
scatter of the M» — Lk relation is still larger than that of the Mm —o relation, while the 
Mm — Mq relation is as tight as the Mm — cr relation. However, the Mm — Mqct 2 relation 
undoubtedly gives the best fit of the experimental data, see Table 1. This means that 
A/qct 2 is a better tracer of Mm than L^, Mq, and a. 
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Table 4 The data are all taken from Hu (2009) [29]. The la errors of L K and a c are 10% and 
5% respectively, whereas that of log Mq are 0.15 dex. 



Galaxy 


Type 


a e (km/s) Log 10 




M. 


(Mq) 


SM, 


.(Mq) 


M G (M ) 


NGC821 


E4 


189 


10.75 


8.5 


x 


10 7 


5.9 


x 


10 7 


4.5 x 


10 10 


NGC1399 


El 


317 


11.35 


1.2 


x 


10 9 


1.2 


x 


10 9 


3.7 x 


10 11 


NGC2974 


E4 


233 


10.95 


1.7 


x 


10 8 


2.1 


X 


10 7 


1.3 x 


10 11 


NGC3377 


E5 


138 


10.18 


1.1 


x 


10 8 


1.1 


X 


10 8 


8.7 x 


10 9 


NGC3379 


El 


201 


10.86 


1.2 


x 


10 8 


1.2 


x 


10 8 


7.4 x 


10 10 


NGC3608 


E2 


178 


10.81 


2.1 


x 


10 8 


1.1 


x 


10 s 


7.1 x 


10 10 


NGC4473 


E5 


192 


10.84 


1.2 


x 


10 8 


3.2 


X 


10 8 


6.3 x 


10 10 


NGC4486 


E0 


298 


11.47 


3.6 


x 


10 9 


1.4 


x 


10 9 


5.9 x 


10 11 


NGC4486A 


E2 


110 


10.26 


1.3 


x 


10 7 


1.5 


x 


10 7 


7.6 x 


10 9 


NGC4552 


E 


252 


11.05 


5.0 


x 


10 8 


6.1 


x 


10 7 


1.1 x 


10 11 


NGC4564 


E6 


162 


10.30 


5.9 


x 


10 7 


8.7 


x 


10 6 


1.8 x 


10 10 


NGC4261 


E2 


309 


11.37 


5.2 


x 


10 8 


1.4 


x 


10 8 


5.0 x 


10 11 


NGC4291 


E2 


242 


10.78 


3.4 


x 


10 8 


9.8 


x 


10 8 


8.3 x 


10 10 


NGC4621 


E5 


211 


10.77 


4.0 


x 


10 8 


4.9 


x 


10 7 


4.4 x 


10 10 


NGC4649 


El 


330 


11.50 


2.0 


X 


10 9 


8.2 


X 


10 8 


5.6 x 


10 11 


NGC4697 


E4 


177 


10.47 


1.7 


X 


10 8 


3.4 


X 


10 7 


3.5 x 


10 10 


NGC5077 


E3 


261 11.26 


7.2 


X 


10 8 


5.1 


X 


10 8 


2.2 x 


10 11 


NGC5576 


E3 


183 


10.50 


1.8 


X 


10 8 


5.2 


X 


10 7 


5.9 x 


10 1() 


NGC5813 


El 


230 


11.06 


7.1 


X 


10 8 


8.6 


X 


10 7 


1.1 x 


10 11 


NGC5845 


E3 


239 


10.54 


2.6 


X 


10 8 


3.7 


X 


10 s 


3.8 x 


10 1() 


NGC5846 


E0 


238 


11.33 


1.1 


X 


10 9 


1.1 


X 


10 8 


3.5 x 


10 11 


NGC6251 


E2 


290 


11.82 


6.2 


X 


10 8 


3.2 


X 


10 8 


6.2 x 


10 11 


NGC7052 


E4 


266 


11.39 


4.0 


X 


10 8 


2.8 


X 


10 8 


3.0 x 


10 11 


IC4296 


E 


322 


12.36 


1.3 


X 


10 9 


2.4 


X 


10 8 


1.9 x 


10 12 


CygnusA 


E 


270 


12.07 


2.9 


X 


10 9 


9.2 


X 


10 8 


2.9 x 


10 12 


IC1459 


E3 


340 


11.54 


2.5 


X 


10 9 


5.1 


X 


10 8 


3.8 x 


10 11 


NGC221 


E2 


75 


3.780 


2.5 


X 


10 6 


6.5 


X 


10 6 


2.8 x 


10 8 


NGC4742 


E4 


90 


10.27 


1.4 


X 


10 7 


7.8 


X 


10 6 


3.5 x 


10 9 



NGC524 


SO 


235 


11.23 


8.3 


X 


10 8 


5.9 x 


10 7 


2.6 x 10 11 


NGC2549 


SO 


145 


10.18 


1.4 


X 


10 7 


5.2 x 


10 7 


1.8 x 10 1() 


NGC3115 


SO 


230 


10.42 


9.3 


X 


10 8 


5.5 x 


10 8 


2.6 x 10 1() 


NGC3245 


SO 


205 


10.64 


2.1 


X 


10 8 


6.6 x 


10 7 


4.1 x 10 1() 


NGC3414 


SO 


205 


10.71 


2.5 


X 


10 8 


3.7 x 


10 7 


5.0 x 10 1() 


NGC3585 


SO 


213 


11.18 


3.4 


X 


10 8 


1.5 x 


10 s 


1.1 x 10 11 


NGC3607 


SO 


229 


11.09 


1.2 


X 


10 8 


6.2 x 


10 7 


1.8 x 10 11 


NGC3998 


SO 


268 


10.85 


2.9 


X 


10 8 


5.8 x 


10 7 


6.9 x 10 1() 


NGC4026 


SO 


180 


10.57 


2.1 


X 


10 8 


6.6 x 


10 7 


3.4 x 10 1() 


NGC4459 


SO 


168 


10.54 


7.1 


X 


10 7 


1.6 x 


10 7 


2.3 x 10 1() 


NGC5128 


SO 


138 


10.33 


5.0 


X 


10 7 


6.1 x 


10 6 


1.3 x 10 1() 


NGC5252 


SO 


190 


11.84 


1.0 


X 


10 9 


1.5 x 


10 9 


1.4 x 10 11 


NGC7457 


SO 


78 


9.64 


3.8 


X 


10 6 


1.7 x 


10 6 


5.4 x 10 9 


P49940 


SO 


288 


11.17 


3.9 


X 


10 9 


5.7 x 


10 8 


7.6 x 10 11 


NGC1023 


SBO 


205 


10.49 


4.4 


X 


10 7 


5.3 x 


10 6 


2.7 x 10 1() 


NGC1316 


SBO 


226 


11.25 


1.6 


X 


10 8 


3.3 x 


10 7 


9.3 x 10 1() 


NGC4596 


SBO 


152 


10.38 


7.8 


X 


10 7 


5.7 x 


10 7 


2.0 x 10 1() 


NGC2787 


SBO 


218 


9.82 


4.1 


X 


10 7 


6.0 x 


10 6 


2.1 x 10 1() 


NGC3384 


SBO 


143 


10.43 


1.7 


X 


10 7 


2.6 x 


10 6 


1.3 x 10 1() 



Circinus 


S 


75 


9.80 


1.1 


X 


10 6 


2.5 


X 


10 6 


2.0 x 10 9 


NGC3393 


S 


184 


10.81 


3.1 


X 


10 7 


2.2 


X 


10 6 


1.0 x 10 11 


NGC224 


S 


160 


10.22 


1.4 


X 


10 8 


9.3 


X 


10 7 


1.9 x 10 1() 


NGC3031 


S 


173 


10.40 


7.9 


X 


10 7 


2.1 


X 


10 7 


1.9 x 10 1() 


NGC4151 


s 


97 


10.27 


3.2 


X 


10 7 


8.2 


X 


10 7 


7.1 x 10 9 


IC2560 


SB 


137 


10.48 


2.9 


X 


10 6 


7.5 


X 


10 6 


2.3 x 10 1() 


Milky Way 


SB 


103 


10.25 


4.1 


X 


10 6 


0.6 


X 


10 6 


1.3 x 10 1() 


NGC1068 


SB 


165 


10.81 


8.3 


X 


10 6 


3.9 


X 


10 6 


3.9 x 10 10 


NGC3079 


SB 


146 


10.46 


2.5 


X 


10 6 


2.5 


X 


10 6 


1.7 x 10 1() 


NGC3227 


SB 


131 


9.89 


2.0 


X 


10 7 


3.8 


X 


10 7 


1.2 x 10 1() 


NGC4258 


SB 


148 


9.93 


3.9 


X 


10 7 


9.1 


X 


10 6 


1.2 x 10 10 
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3 Age— temperature diagram for galaxies 

The main finding that emerges from the analysis of the data is that the values of 
the intrinsic scatter and of the \ 2 °f the relation between the mass of SMBH located 
in the center of the host galaxies and the kinetic energy of random motions of the 
corresponding bulges are better than those of the other relations. This statement is 
true for each of the three samples considered in this work, which have been extracted 
from independent catalogues. Besides having a support from numerical models [17,18], 
the M, — Mqo 1 law is also the only one with a quite clear physical explanation. In 
fact, the mass of the central BH, just like the entropy, can only increase with time or 
at most remain the same but never decrease. So, M» is connected with the age of the 
galaxy. On the other side, the kinetic energy of the stellar bulges directly determines 
the temperature of the galactic system [26]. Hence, it is possible to reinterpret our 
relation as an age-temperature diagram for galaxies. 

There is also an evident indication that the galaxies occupy quite well-defined 
regions of the M» — (Mg<t 2 )/c 2 plane, accordingly to their morphological type. In fact, 
they are distributed along the best-fitting line, in such a way that the ellipticals are 
in the upper part of the diagram, the lenticulars in the middle, and the spirals in the 
lower part, see Fig. 1(c), Fig. 2(d), and Fig. 3(d). The distribution of the galaxy types 
changes a little bit moving from one sample to another depending on the differences 
of the morphological classification of some galaxies given by the various authors of 
the catalogues. This is essentially caused by the fact that whenever the lenticulars are 
inclined face-on it is often difficult to distinguish between them and ellipticals. On the 
other hand, many of the flattest ellipticals may actually be misclassified lenticulars [57]. 
Actually, most of the lenticulars are spread and mixed with the low/intermediate- 
mass ellipticals, supporting the emerging idea that their nature is different from the 
standard one that sees them just as a transition type between ellipticals and spirals [58]. 
The barred lenticulars are generally located in a lower zone than that occupied by 
lenticulars. The few galaxies of types E4 and E6 appear only in the middle region, this 
because the mass of their SMBH is roughly an order of magnitude lower than than 
those in host galaxies of types EO and El. Finally, spirals and barred spirals are located 
in the lower zone of the FM diagram, even if their distribution above and under the 
best-fitting line is not clear since it changes dramatically depending on the sample 
considered. 

Together with the spirals, we find some dwarf elliptical galaxies (dEs) in the lower 
part of the diagram, that supports the hypothesis that the dEs could have suffered a 
morphology changing due to a process called "galaxy harassment" [59]. To be more 
precise, the dEs could be the remnants of low-mass spiral galaxies that obtained a 
rounder shape through the action of repeated gravitational interactions with giant 
galaxies within a cluster. If this theory is exact, than it is possible to explain the 
presence of the dEs in the zone of the FM diagram where most of the late type galaxies 
are confined. On the other hand, the idea that the dEs may be primordial objects, rather 
than young galaxies, is in part supported by numerical simulations, which predict a 
very low SMBH mass limit for the ellipticals, smaller with respect to observations, as 
is discussed in the next section. This hypothesis could be correct if new ellipticals will 
be found in the middle region of the diagram. 

For early type (E + SO) galaxies, the ratio of the galaxy mass to SMBH mass 
is essentially independent of galaxy flattening, i.e. the fraction of the total mass of 
early-type galaxies that is in the form of a central BH is independent of the angular 



12 



momentum of the host galaxy. In fact, using the data of the galaxy sample of Curir et 
al. [60], the correlation coefficient of the relation between the ratio (M./Mq) and the 
specific angular momentum is very low, r = 0.54. This agrees with the fact that the 
kinetic energy of stellar random motions of the early type galaxies is higher than that 
of the late type galaxies (see the differences between the values of a in Tables 2, 3, 
and 4). However, it is puzzling that the lenticular galaxies NGC4342 and NGC4350 
present a high values of the SMBH mass to bulge mass ratio, see Fig. 1(c). Their 
distance from the best-fitting line is quite suspicious and maybe is an indication of a 
peculiar evolutionary story. 

4 Numerical models 

We compare the observed age-temperature relation described in the previous sections 
with the predictions of two hierarchical models of structure formation, one by De 
Lucia & Blaizot 2007 [61] (the 'Max Planck Institut fur Astrophysik (MPA) model') 
and the other by Bower et al. 2006 [62] (the 'Durham model'), that simulate the 
cosmological co-evolution of dark matter haloes, subhaloes, galaxies and SMBH in 
the ylCDM-cosmology framework. Both models have been implemented on top of a 
high-resolution cosmological N-body simulation, the Millennium Run [63], that follows 
the evolution of N = 2160 3 dark matter particles of mass 8.6 x 1O 8 /i _1 M0, within 
a co-moving box of size 500/i _1 Mpc on a side and scale resolution of 5 h 1 kpc , 
from 2 = 127 to the present. The cosmological parameters adopted are O m — 0.25, 
J?b = 0.045, h = 0.73, O a = 0.75, n = 1, and a s = 0.9 [64]. The output of both 
models, dubbed DeLucia2006a and Bower2006a, respectively, are publicly available at 
http://www.mpa-garching.mpg.de/millennium [65]. 

The high mass and spatial resolution of the Millennium Simulation allows one to 
track the motion of dark matter substructures inside massive haloes, making it possible 
to construct merging history trees of all the dark matter haloes and subhaloes inside 
the simulation box [63]. The dynamical evolution of all satellite galaxies is followed 
until tidal truncation and stripping disrupt their host dark matter subhaloes. Then a 
residual survival time is estimated by computing the dynamical friction formula. To 
populate the dark matter subhaloes with galaxies and black holes, both models adopt 
a set of equations to describe the radiative cooling of gas, the star formation, the 
metal enrichment and supernovae feedback, the growth and feedback of SMBH, the 
UV background reionization, and the effects of galaxy mergers. 

In the MPA model, the BH mass accretion is triggered by two different phenomena: 
i) the merger between gas-rich galaxies (quasar mode) and ii) the cooling flow at the 
centers of X-ray emitting atmospheres in galaxy groups and clusters (radio mode). In 
the Durham model, in addition to these two mechanisms, the BH mass accretion is also 
triggered by the disk instability. The most important consequence of the BH growth 
in the evolution of galaxies is that the energy feedback associated with the radio mode 
accretion is assumed to be able to reduce or even stop the cooling flow in the halo 
centers. 

A full description of the MPA and Durham models and a comparison between 
their main predictions with observations can be found in [66,61,18,67,68] and in [62, 
69], respectively. 

In both models the morphology of a galaxy is determined by the bulge-to-disc ratio 
of its absolute, rest frame B-band luminosity. Specifically, a galaxy is classified as early 
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Fig. 4 The local M. — Mqct 2 relation predicted by the two cosmological models considered. 
Red, green and blue dots refer to early type, SO and late type galaxies, respectively. The black 
lines show the linear best fit relations to the observational dataset, reported in Table 1: Graham 
2008 sample [27] - solid line, Giiltekin et al. 2009 sample [28] - dashed line and Hu 2009 [29] 
sample - dotted lines. Only 1/100 of the total galaxy sample is shown in the figure. 
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Fig. 5 The local M, — Mqct 2 relation predicted by the MPA model in a compact form, to 
directly compare with Fig. 1(c), 2(c) and 3(c). The symbols are as is Fig. 4. 



type if AM < 0.4, as lenticular (SO) if 0.4 < AM < 1.56 and as late type if AM > 1.56, 
where AM — Afbulgo — -^total and M|-, u ig e and M tota i are the B-band magnitude of 
the bulge and of the whole galaxy, respectively [61]. 

In Fig. 4 we show the local (z = 0) M» — Mq(j 2 relation predicted by the two models 
described above. As explicitly indicated by the labels, red, green and blue dots refer to 
early type, SO and late type galaxies, respectively. The bulge velocity dispersion, a, is 
derived from the V c — a relation of Baes et al. 2003 [70] , as in [18] . To closely mimic the 
observed catalogues, we have selected the model galaxies in these three observed mass 
ranges: 3.7 x 10 10 < M G /M < 1.6 x 10 12 for early type, 10 10 < M G /Mq < 2.4 x 10 11 
for SO, and 1.7 x 10 9 < M G /M© < 1.3 x 10 11 for late type galaxies. Besides, for 
clarity reasons, we have shown only 1/100 of the total galaxy sample. The black lines 
show the linear best fit relations to the observational dataset, reported in Table 1: 
Graham 2008 sample [27] - solid line, Giiltekin et al. 2009 sample [28] - dashed line 
and Hu 2009 [29] sample - dotted lines. As can be seen in Fig. 4, the predictions of 
both the MPA and Durham model are in quite good agreement with the observed 
scaling relations. Also, the upper BH mass limits of the three groups considered (early 
type, SO and late type) match the observed ones, while the lower BH mass limits 
predicted are smaller with respect to observations. More accurate observational data 
are necessary to understand if this is a problem of the models themselves or simply 
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depends on observational selection effects due to the fact that the number of galaxies 
with a realiable measured SMBH mass is still low. 

Finally, to directly compare the MPA model predictions with the observational 
data reported in Fig. 1(c), 2(d) and 3(d), we have shown in Fig. 5 the M. — Mqo 1 
relation in a more compact form. 

5 Summary 

We have tested the goodness of the M. — Mqo 1 relation as a predictor of the SMBH 
mass in the center of galaxies on three different samples based on the galaxy catalogues 
of Graham [27], Giiltekin et al. [28], and Hu [29]. Our analysis shows that this relation 
is the one with the lowest scatter when compared with other relations like M. — a, 
M» — Mq, and M, — Lq. This is evident if we look into the figures 1, 2, and 3, and the 
Table 1, where the values of the intrinsic scatter and of the x 2 are reported, providing 
a quick comparison among the various relations. The galaxies arrange themselves on 
the M. — Mq<t 2 plane following the best-fitting line and in accordance with their 
morphological type. In particular, the late types are clearly separated from the early 
types; the lenticulars are more mixed with the ellipticals; the barred lenticulars are 
more mixed with the spirals; the dwarf ellipticals occupy the same region of the spirals 
possibly revealing themselves as remnants of low-mass spirals, see Fig. 1(c), 2(d) 
and 3(d). From this point of view, the M. — Mqcj 2 relation can be revised as an age- 
temperature diagram for the galaxies, where the age is deducted from the mass of the 
central SMBHs, while the temperature from the kinetic energy of the random motions 
of the stars belonging to the spheroidal components [26]. 

The soundness of the FM diagram has also been investigated in the ylCDM cosmol- 
ogy using two galaxy formation models based on the Millennium Simulation, one by 
Bower et al. (the Durham model) [62] and the other by De Lucia & Blaizot (the MPA 
model) [61]. Both the MPA and the Durham model reproduce quite well the real data, 
taking into account the uncertainties of the presentday observational measurements. 
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